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Abstract 

We document the discovery of two generating functions for £(2n + 2), analogous to 
earlier work for £(2n + 1) and £(4n + 3), initiated by Koecher and pursued further by 
Borwein, Bradley and others. 

1 Introduction 

Stimulated by recent work in the arena of Apery- like sums 13 Ej we decided to method- 
ically look for series acceleration formulas for the Riemann zeta function involving central 
binomial coefficients in the denominators. Using the PSLQ integer relation algorithm, as 
described below, we uncovered several new results. In particular, we document the discov- 
ery of two generating functions for C(2n + 2), analogous to earlier work for ((2n + 1) and 
£(4n + 3), initiated by Koecher and pursued further by Borwein, Bradley and others. As a 
conclusion to a very satisfactory experiment, we have been able to use the Wilf-Zeilberger 
technique to prove our results. 

An integer relation detection algorithm accepts an n-long vector x of real numbers 
and a bound A as input, and either outputs an n-long vector a of integers such that the 
dot product a\X\ + • • • + a n x n = to within the available numerical precision, or else 
establishes that no such vector of integers of length less than A exists. Here the length is 
the Euclidean norm (a\ + af, + • • • + a 2 ) 1 / 2 , derived from the usual Euclidean metric on 
R n . Helaman Ferguson's PSLQ |131 12] is currently the most widely used integer relation 
detection algorithm pp. 230-235], although variants of the so-called LLL algorithm |16j 
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are also commonly employed. Such algorithms underlie the "Recognize" and "identify" 
commands in the respective computer algebra packages Mathematica and Maple. They 
also play a fundamental role in the investigations we discuss here. 

This origins of this work lay in the existence of infinite series formulas involving central 
binomial coefficients in the denominators for the constants C(2), C(3), and C(4). These 
formulas, as well the role of the formula for ("(3) in Apery's proof of its irrationality, have 
prompted considerable effort during the past few decades to extend these results to larger 
integer arguments. The formulas in question are 

C(2) = 3 g^(f)' W 
CP) - 2^V^P < 2 > 



k 

00 



cw - % E^- (3) 

Identity Q has been known since the 19th century — it relates to arcsin 2 (x) — while P|) 
was variously discovered in the last century and © was noted by Comtet p. 89], 
see [SJ E3- Indeed, in j^j a coherent proof of all three was provided in the course of a 
more general study of such central binomial series and so-called multi-Clausen sums. 
These results led many to conjecture that the constant Q5 defined by the ratio 

/ 00 (—■nk+i 

a - c<5) /£w 

is rational, or at least algebraic. However, integer relation computations using PSLQ and 
10,000-digit precision have established that if Q5 is a zero of a polynomial of degree at 
most 25 with integer coefficients, then the Euclidean norm of the vector of coefficients 
exceeds 1.24 x 10 383 . Similar computations for ("(5) have yielded a bound of 1.98 x 10 380 . 

These computations lend credence to the belief that Q5 and £(5) are transcendental. If 
algebraic, they almost certainly satisfy no simple polynomial of low degree. In particular, 
if there exist relatively prime integers p and q such that 

p CO { _ l)k+ l 



then p and q must be astronomically large. Moreover, a study of polylogarithmic ladders 
in the golden ratio produced ^ 

2 ^ 5 )-E^?7^ = 2 Li 5^)-2 Li4 ^ )l0g ^ + C(3)l0g2 ^-3 C(2)l0g3p -^ l0g5 ^ 

k=l \ k J 

(4) 
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where p = (3 — \/5)/2 and hi n (z) = Ylk=i z V^ n ls the polylogarithm of order n. Since the 
terms on the right hand side of (j3J) are almost certainly algebraically independent _TT] , we 
see how unlikely it is that Q5 is rational. Although the irrationality of £(5) has not yet 
been confirmed, it is known that one of C(5), C(7), C(9)> C(H) is irrational |2U| . 

Given the negative result from PSLQ computations for Q5, the authors of [7j systemat- 
ically investigated the possibility of a multi-term identity of this general form for £(2n+l). 
The following were recovered early |Hj in experimental searches using computer-based 
integer relation tools: 

c(5) = 2 Sw^w^' (5) 

5 (-!)*+* 25 (-!)*+! *Z* 1 

c(9) = 4g^^"ig^frg^ +5 S^(frS^ 

45^ (-l) fc+1 ^l 25 ^ (-l) fc+1 ^ 1 ^ 1 
fc=l K \k) j=i J k=l K \k) j=l J j=l J 



5A(^ 25A(-ir^l 

1 j ^fe ,t^ 4 



75 ^ (-l) fc+1 1 125 ^ (-l) fc+1 ^1^1 



The general formula 



1 1 ^ (-l) fc+1 5A: 2 - x 2 ^r 1 / x 2 \ 

^ k(k 2 -x 2 ) = 2^ fc3(2^ fc 2 -x 2 11 l 1 "^ 2 ) ' <); ' 



fc 7 m=l 



was obtained by Koecher ^3] following techniques of Knopp and Schur. It gives © as 
its first term and © as its second term but more complicated expressions for C(7), C(9) 
and C(H) than (0) and (JHJ). The corresponding result that gives ©, © and © for 
its first three terms was worked out by Borwein and Bradley 0. 

Using bootstrapping and an application of the "Pade" function (which in both Math- 
ematica and Maple produces Pade approximations to a rational function satisfied by a 
truncated power series) produced the following remarkable and unanticipated result [7j: 

2-F(l-^A 4 ) 2^ i k*[™)(l- x ±/k*)li 1 \l-xVrn* J' 1 ] 
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The equivalent hyper geometric formulation of (|1(J[) is 



5^4 



2, 1 + x, 1 — x, 1 + ix, 1 — ix 
2 + x, 2 — x, 2 + ix, 2 — ix 

5 



6-^5 



2, 2, 1 + x + ix, 1 + x — ix, 1 — x + ix, 1 — x — ix 
3/2,2 + x,2 - x,2 + ix,2 - ix 



The identity (fTU)) generates ©, (jHJ) and © above, and more generally gives a formula 
for £(4n + 3), which for n > 1 contains fewer summations than the corresponding formula 
generated by @. The task of proving (|TU)) was reduced in |7j to that of establishing any 
one of a number of equivalent finite combinatorial identities. One of these latter identities 
is 



£^n(« 4 +' 4 >/ri< 

k=l i=l 1 i=l 

i^k 



; 4 - i 4 l 



2n 
n 



(11) 



This was proved in PQ, so (|1U|) is an established theorem. It is now known to be the x = 
case of the even more general formula 



E 

k=l 



k 



fc 4 — x 2 k 2 



1 ~ (-1 

9 7./: 



,fc+l 



5k 2 



c 2 /c 2 



n 

m=l 



x 2 ) 2 + 4y 4 



m 4 — x 2 m 2 



, (12) 



in which setting y = recovers @. The bivariate generating function identity 1)12(1 
was conjectured by Henri Cohen and proved by Bradley |10| . It was subsequently and 
independently proved by Rivoal [T8| . 

Following an analogous — but more deliberate — experimental-based procedure, as de- 
tailed below, we provide a similar general formula for C(2n. + 2) that is pleasingly parallel 
to (UUJ). It is: 



Theorem 1 Let x be a complex number not equal to a non-zero integer. Then 



E 



1 



k 2 



1 



fc-i 

n 



1 — 4x 2 /m 2 

k=i - - k=i k 2 ( 2 k k ) (1 - x 2 /fc 2 ) V 1 " x 2 /m 2 

Note that the left hand side of (|13|) is trivially equal to 

f>n + 2)z 2 " = 1"™^) 

n=0 



(13) 



2x 2 



(14) 



Thus, l)13|) generates an Apery- like formulae for <^(2n) for every positive integer n. 

In Section 2 we shall outline the discovery path, and then in Section 3 we prove (fT3|l - 
or rather the equivalent finite form 



3-^2 



3n, n + 1, — n 
2n + l,n + 1/2 



( 3 ;) 



< n G Z. 



(15) 



In Section 4 we provide another generating function for which the leading term or 
"seed" is Comtet's formula (J2J) for C(4), while the prior generating functions have seeds 
© and ©. 

The paper concludes with some remarks concerning our lack of success in obtaining 
formulas analogous to (|1U|1 and which would generate the simplest known Apery-like 
formulae for £(4n + 2) and ((An + 1), respectively. In this light we record 

^ 1 2tpv/3 + 9 nr , 

Z^72feV = — 2l — ' ( } 



k=l \k. 



which are perhaps more appropriate seeds in these cases, see pp. 384- 



2 Discovering Theorem [T] 

As indicated, we have applied a more disciplined experimental approach to produce an 
analogous generating function for Q(2n-\-2). We describe this process of discovery in some 
detail here, as the general technique appears to be quite fruitful and may well yield results 
in other settings. 

We first conjectured that ((2n + 2) is a rational combination of terms of the form 

oo ^ N k—1 

a(2r-[2 ai ,...,2a N }) := Et^IIE^' ( 18 ) 

k=l \ k I i=l rij=l i 

where r + a i = n + 1, and the o« are listed in nonincreasing order (note that the 
right-hand-side value is independent of the order of the Oj). This dramatically reduces the 
size of the search space, while in addition the sums (|18|) are relatively easy to compute. 
One can then write 

oo oo n+1 

EC(2n + 2)x 2 " = EE E «(vr)a(2r;27r)x 2 ™, (19) 

n=0 n=0 r=l 7ren(n+l-r) 

where II(m) denotes the set of all additive partitions of m if m > 0, 11(0) is the singleton 
set whose sole element is the null partition [], and the coefficients a(ir) are complex 
numbers. In principle a(n) in ()19|) could depend not only on the partition it but also on 
n. However, since the first few coefficients appeared to be independent of n, we found it 
convenient to assume that the generating function could be expressed in the form given 
above. 

For positive integer k and partition it = (aa2, ■ ■ ■ ,ajv) of the positive integer m, let 

N k-1 

n e i^T' 

i=l n,=l i 



Sfc(7r) := 



5 



Then 



C r(2r;27r) = >T 



gfc(7[) 



and from ((TU|) . we deduce that 

oo oo n+1 

^C(2n + 2)x 2 ™ = EE E «(^(2r;2vr)^ 

n=0 n=0 r=l 7r€ll(n+l-r) 



EyivE V E E -w^w^ 

fc=l IfeJ r=l n=r-l 7ren(n+l-r) 

OO -. oo 

E n(fc2 x2) E* 2m E 

H Ui^ X > m=0 7rGn(m) 

E 72^ iT^Cx) (20) 

k=l {k)( k - x ) 



where 



flfc(x) := J]^ £ a(7r)a fc (vr), (21) 

m=0 7rGn(m) 

whose closed form is yet to be determined. Our strategy, as in the case of (fTU)l 8j, was to 
compute Pk{x) explicitly for a few small values of k in a hope that these would suggest a 
closed form for general k. 

Some examples we produced are shown below. At each step we "bootstrapped" by 
assuming that the first few coefficients of the current result are the coefficients of the 
previous result. Then we found the remaining coefficients (which are in each case unique) 
by means of integer relation computations. 

In particular, we computed high-precision (200-digit) numerical values of the assumed 
terms and the left-hand-side zeta value, and then applied PSLQ to find the rational 
coefficients. In each case we "hard-wired" the first few coefficients to agree with the 
coefficients of the preceding formula. Note below that in the sigma notation, the first 
few coefficients of each expression are simply the previous step's terms, where the first 
argument of a (corresponding to r) has been increased by two. 

These initial terms (with coefficients in bold) are then followed by terms with the other 
partitions as arguments, with all terms ordered lexicographically by partition (shorter 
partitions are listed before longer partitions, and, within a partition of a given length, 
larger entries are listed before smaller entries in the first position where they differ; the 
integers in brackets are nonincreasing) : 



6 



OO -. 

C( 2 ) = 3 E72^ = 3 ^ 2 >[°])> ( 22 ) 

fc=i U ) K 

oo i oo v-^fc — 1 -—2 

C(4) = 3£^-9£^— = M4,[0])-M2,[2]) (23) 



ESr 2 45 ~ eHj- 4 



27^^E?=ii- 2 



+ 2 EE A2(i\ k 2 ' ^ 24 ^ 

fc=l j=l J \k) K 

45 27 

= 3a(6, []) - 9<r(4, [2]) - T a(2, [4]) + y a(2, [2, 2]) (25) 

C(8) = 3<x(8, []) - 9a(6, [2]) - ^a(4, [4]) + ^a(4, [2, 2]) - 63a(2, [6]) 

135 2T 

+— <r(2,[4,2])- y <7(2,[2,2,2]) (26) 

45 27 

C(10) = 3(7(10, []) - 9a(8, [2]) - -^(6, [4]) + —a(6, [2, 2]) - 63a(4, [6]) 

iqi; 97 7fi5 

+^a(4, [4, 2]) - ^-a(4, [2, 2, 2]) - ^a(2, [8]) + 189a(2, [6, 2]) 

+^a(2, [4, 4]) - ?fa(2, [4, 2,2]) + fa(2, [2, 2, 2,2]). (27) 

Next from the above results, one can immediately read that «([]) = 3, a([l]) = —9, 
a([2]) = —45/2, a([l,l]) = 27/2, and so forth. Table 1 presents the values of a that we 
obtained in this manner. 

Using these values, we then calculated series approximations to the functions P&(x), 
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Partition 


Alpha 


Partition 


Alpha 


Partition 


Alpha 


[empty] 


3/1 


1 


-9/1 


2 


-45/2 


1,1 


27/2 


3 


-63/1 


2,1 


135/2 


1,1,1 


-27/2 


4 


-765/4 


3,1 


189/1 


2,2 


675/8 


2,1,1 


-405/4 


1,1,1,1 


81/8 


5 


-3069/5 


4,1 


2295/4 


3,2 


945/2 


3,1,1 


-567/2 


2,2,1 


-2025/8 


2,1,1,1 


405/4 


1,1,1,1,1 


-243/40 


6 


-4095/2 


5,1 


9207/5 


4,2 


11475/8 


4,1,1 


-6885/8 


3,3 


1323/2 


3,2,1 


-2835/2 


3,1,1,1 


567/2 


2,2,2 


-3375/16 


2,2,1,1 


6075/16 


2,1,1,1,1 


-1215/16 


1,1,1,1,1,1 


243/80 


7 


-49149/7 


6,1 


49140/8 


5,2 


36828/8 


5,1,1 


-27621/10 


4,3 


32130/8 


4,2,1 


-34425/8 


4,1,1,1 


6885/8 


3,3,1 


-15876/8 


3,2,2 


-14175/8 


3,2,1,1 


17010/8 


3,1,1,1,1 


-1701/8 


2,2,2,1 


10125/16 


2,2,1,1,1 


-6075/16 


2,1,1,1,1,1 


729/16 


1,1,1,1,1,1,1 


-729/560 


8 


-1376235/56 


7,1 


1179576/56 


6,2 


859950/56 


6,1,1 


-515970/56 


5,3 


902286/70 


5,2,1 


-773388/56 


5,1,1,1 


193347/70 


4,4 


390150/64 


4,3,1 


-674730/56 


4,2,2 


-344250/64 


4,2,1,1 


413100/64 


4,1,1,1,1 


-41310/64 


3,3,2 


-277830/56 


3,3,1,1 


166698/56 


3,2,2,1 


297675/56 


3,2,1,1,1 


-119070/56 


3,1,1,1,1,1 


10206/80 


2,2,2,2 


50625/128 


2,2,2,1,1 


-60750/64 


2,2,1,1,1,1 


18225/64 


2,1,1,1,1,1,1 


-1458/64 


1,1,1,1,1,1,1,1 


2187/4480 











Table 1: Alpha coefficients found by PSLQ computations 
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by using formula (|21|). We obtained: 



P 3 (x) 



45 2 
45 



45 
16' 



-.x 



45 
64 : 



45 
256' 



-x 



45 

1024 : 



.10 



45 
4096 ^ 



,12 



45 



.14 



16384 



.16 



65536 



49 9 119 4 3311 4 38759 

3 x 2 H x 4 H H 

4 144 5184 186624 



384671 
6718464' 



+ 



3605399 
241864704' 



,.io 



33022031 
8707129344' 



.12 



205 



Ps{x) ~ 3-— r 



7115 



16 



2304 



x* + 



207395 
331776 



x u + 



299492039 

313456656384 : 
4160315 



x 8 + 



.14 



74142995 



-x 



10 



+ 



1254489515 
990677827584 



x 12 + 



x 14 + 



5269 2 6640139 



400 



1440000 



20685646595 

142657607172096 

4 1635326891 „ 

x + x 



5184000000 



47775744 6879707136 
336494674715 

20542695432781824 ; 
5944880821 8 

18662400000000 X 



,16 



212874252291349 
67184640000000000 ' 



.10 



141436384956907381 
241864704000000000000 ' 



,12 



70524260274859115989 
870712934400000000000000 ; 



.14 



31533457168819214655541 
3134566563840000000000000000 : 



.16 



P 7 (x) 



5369 2 8210839 

~ 400 X + 1440000 ^ 
278500311775049 

67184640000000000 ' 



199644809 
5184000000 : 



680040118121 
18662400000000 : 



,10 



84136715217872681 
241864704000000000000 ' 



,12 



22363377813883431689 
870712934400000000000000 ' 



,14 



5560090840263911428841 
3134566563840000000000000000 ' 



.16 



With these approximations in hand, we were then in a position to attempt to determine 
closed- form expressions for Pk(x). This can be done by using either "Pade" function in 
either Mathematica or Maple. We obtained the following: 



Pi(x 
P 2 (x 

P 3 (x 

P 4 (x 

P 5 (x 

P 6 (x 

Prix 



3 

3(4x 2 



1) 



(x 2 - 1) 
12(4x 2 - 1) 

(x 2 - 4) 
12(4x 2 - l)(4x 2 



9) 



48(4x 2 



4)(x 2 - 9) 
- l)(4x 2 - 9) 



(x 2 -9)(x 2 - 16) 
48(4x 2 - l)(4x 2 -9)(4x 2 



25) 



(x 2 - 9)(x 2 - 16) (x 2 - 25) 
192(4x 2 - l)(4x 2 - 9)(4x 2 - 25) 
(x 2 - 16) (x 2 - 25) (x 2 - 36) 
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These results immediately suggest that the general form of a generating function iden- 
tity is: 



k - 1 < „2 _2 



£ C (2» + 2)x- I s^-^-L—n^L, {28) 



I (k 2 — x 



n=0 k=l \k Jy I m=l 

which is equivalent to (|13(1 . 

We next confirmed this result in several ways: 

1. We symbolically computed the power series coefficients of the LHS and the RHS of 
(|28|) . and have verified that they agree up to the term with x 100 . 

2. We verified that £(1/6), where Z(x) is the RHS of (j2SJ), agrees with 18 - 3^/3^, 
computed using (|14|). to over 2,500 digit precision; likewise for £(1/2) = 2, £(1/3) = 
9/2 - 3vr/(2V3), £(1/4) = 8 - 2vr and £(1/^2) = 1 - vr/v 7 ^ • cot(vr/ v / 2). 

3. We then affirmed that the formula (|28[) gives the same numerical value as 1)14(1 for 
the 100 pseudorandom values {mir}, for 1 < m < 100, where {•} denotes fractional 
part. 

Thus, we were certain that (|13|) was correct and it remained only to find a proof of 
Theorem 1. 



3 Proof of Theorem U 

By partial fractions, 

1 ^r 1 1 - 4z 2 /m 2 _ A cn(fc) 

' m=l ' n=l ' 

where 



fe-i 



Cn(k) = \[{l-An 2 /m 2 ) \\{l-n 2 / 



m 2 ) 



m=l 1 m=l 



if 1 < n < k, and c n (/c) = if n > A; or if k > 2n + 1. It follows that 

1 ^ 1 _ 4z 2/ m 2 - J * Cn{k) 



v - 1 TT 1 ~ 4z z /m z v - 1 x - 

^F^a-^iA 1 — z 2 /m 2 " ^^J^I-^/n: 



^2/^2 

oo 2n ,,s 

El C n {K) 
1 - z 2 /ri 2 ^ 

n=l ' fc=n * I k ) 



The interchange of summation order is justified by absolute convergence. To prove (|13|) . 
it obviously suffices to show that 

2n 2ri fc— 1 / k 

k=n \ k ) k=n V k J m=l ' m=l 
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for each positive integer n. But, 



6n 



(2n)! 



2 2n 1 fc— 1 / fc 

J Ein^ 2 -- 2 )/ n (- 2 -- 2 ) 



k=n 
In! 



2fc\ 
. fc J 



m=l 



m=n+l 



(3n)! n 
(2ra)!(2ra) 



3n, n + 1, — n 



3i?2k 2n + l,ra + l/2 



Thus, we have reduced the problem of proving (|13j) to that of establishing the finite 
identity 



. = (3n)!n! p/^ 3n,n + l,-n 
1 ' " (2n)!(2n)! 3 2 V 2n + l,n + l/2 



1, 



(29) 



But Maple readily simplifies T(n + 1)/T(n) = 1, and since T(0) = 1, the identity l(2^|) 
and hence ()15j) and ()13j) are established. If a certificate is desired, we can employ the 
Wilf-Zeilberger algorithm. In Maple 9.5 we set 



r :- 



(2n\ 



{3n\ 



f:- 



(3n + k - 1)! (n + k)\ (-n - 1 + A;)! (2n)! (n - 1/2)! (1/4)* 



(3n - 1)! n\ {-n - 1)! (2n + fc)! (n - 1/2 + fc)! fc! 
Maple interprets the latter in terms of the Pochammer symbol 



(30) 



as 



/ 



(«)fc == II( a+ ^- 1 ) 



(3n) fc (n + l) fc (-n)fc (l/4) fc 
(2n + l) fc (n + l/2) fc ' fc! 



so despite the appearance of (|3*0|l the issue of factorials at negative integers does not arise 
for non-negative integers fc and n. Now execute: 

> with(SumTools [Hypergeometric] ) : 

> WZMethod(f ,r,n,k, 'certify') : certify; 

which returns the certificate 

/ 2 \ 
\11 n +l+6n+k+5kn/k 



3 (n - k + 1) (2 n + k + 1) n 
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This proves that summing f(n,k) over k produces r(n), as asserted. 

Indeed, the (suppressed) output of 'WZMethod' is the WZ-pair (F, G) such that 

F(n+l,k)-F(n,k) = G(n,k + 1) - G(n,k), 

where F(n, k) := f(n, k)/r[n) for r(n) ^ and is f(n, k) otherwise. Sum both sides over 
k € Z and use the fact that by construction G(n, k) — > as k — ► ±oo. The certificate is 

G(n,fc) 

K{n,k) := — — — . 
v ' ; F(n, k) 

QED 



4 An Identity for £(2n + 4) 



We compare (f 13|) to a result due to Leshchiner ^7] which is stated incorrectly in [Q, and 
which, as the authors say, has a different flavor: for complex x not an integer, 

if- i 3fc2 + * 2 n7i x2 ) = y t- 1 )"" 1 - l 1 f3n 

2^ k 2(2k\ /c 2 -x 2 m m 2 / ^n 2 -x 2 2xsin(7rx) 2x 2 ' 1 j 

fc=l V fc / m=l n=l v ' 

Using the methods of the previous section — but using a basis of sums over simplices not 
hypercubes — we have likewise now obtained for complex x not an integer, 

1 tt L _ x^_\ k 1 3cos(ttx/3) 

k 2 ( 2 j?)(k 2 - x 2 ) J^Jj V m V 4x 3 sin(7rx) x 4 4x 4 ' 



To see this, let 



i 1/2 \ / 1,1 + a:, l + x,l-x, l-x 
W{X) = (l^) 5Fi { 2,3/2,2 + x,2-x 



denote the left hand side of (|3*2*|) . and let 



V(x) := - 3^2 



i\ / i,i+x,i- x 

2J 32 \ 2,3/2 



i \ °° 1 fe— l / 9 

l \ v A l T — r ( X Z 



fc=l V fc / m= 

2 l (2arcsin(l/2)) 2ri _ l-cos(vrx/3) 



fe=l V fe 7 m=l 

00 /'o /1 /oY\2n 



(2n)! x 2 

n=l 



Expanding Leshchiner's series (|31|) now gives 

3 V(x) + 2x 2 W{x) 



2 2xsin(-7rx) 2x 2 

Solving for W(x) gives (|32|) as claimed. 
To recapitulate, we have 
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Theorem 2 Let x be a complex number, not an integer. Then 

El ^-r^ / x 2 \ nx csc(irx) + 3 cos(7rx/3) — 4 



A': 



If < |x| < 1, toen f/ie Maclaurin series for the left hand side of \'J'J\) is equal to 
(-l) n {3 1 " 2 " - 2B 2n (2 2 "" 1 - 1) } ir 2n 



E 



4(2n)! 



17„.,. 313 , 2 23147 w . 4 1047709 >7 , fi 
= 36 C <"> + 648 C(6 > 1 + 46656 C < 8) 1 + 2099520 C(10) * + ' ' ' • 

where the rational coefficients i?2n refer to i/ie even indexed Bernoulli numbers generated 
by 

^ (2x) 2n 
xcoth(x) = ^-—-B^. 

n=0 V 

Note that the constant term recaptures (j2J) as desired — as taking the limit on the 
right side of (|33|) confirms. Correspondingly, the constant term in (|31|) yields Iflj). The 
coefficient of x 2 is 

313 °° 1 °° 1 k — ^ 1 

This all suggests that there should be a unifying formula for our two identities — (|13j) 
and (|3*3*)) — as there is for the odd cases, see 1)12(1 . 

5 Conclusion 

We believe that this general experimental procedure will ultimately yield results for many 
other classes of arguments, such as for £(4n + m), m = 0, 1, but our current experimental 
results are negative. 

1. Considering £(4n + 1), for n = 2 the simplest evaluation we know is ([7)l. This is one 
term shorter than that given by Rivoal ^Sl , which comes from taking the coefficient 
of x 2 y 4 in (PI) . 

2. For C(2n + 4) (and £(4n)) starting with (j5J) which we recall: 



, , - 36-1 ^-r 1 
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the identity for £(6) most susceptible to bootstrapping is 



C(6) 



36-8 
163 



oo 1 „ oo fe— 1 1 

El O ^ — A 1 ^ — A 1 



For £(8), £(10) we have enticingly found: 

k— 1 -, „ oo ., fc— 1 



C(8) 

C(io) 



36-64 
1373 

36-512 
11143 



y> 1 1 y4 1 1 1 

h k8 ( 2 k k ) + 4 h *Tk) h j4 + ~ 2 k * 2 ( 2 k k ) h j \ 

oo 1 „ oo ^ fc — 1 1 Q oo 1 fc — 1 _. 

A- 1 J-l 

~4 



fc=l V fc / j=l J k- 



9 ^ 1 ^1 27 ^ 1 ^1^1 



fc: =1 j = l J ' 1=1 ^ 



(34) 



(35) 



(36) 



But this pattern is not fruitful; the pattern stops at n = 10. 
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